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The  Rayleigh-Taylor  Instability  in  Ablatively  Accelerated 
Targets  with  1,  1/2  and  1/4  #im  Laser  Light 

I.  INTRODUCTION 

The  Rayleigh-Taylor1  (RT)  hydrodynamic  instability  can  occur  whenever  a  dense  fluid 
is  supported  by  a  less  dense  fluid  in  a  uniform  gravitational  Held;  or  equivalently,  a  dense 
fluid  is  accelerated  by  a  less  dense  fluid.  Small  perturbations  at  the  interface  between  the 
two  fluids  grow  exponentially  in  time  as  the  denser  fluid  interchanges  positions  with  the 
less  dense  fluid.  This  interchange  process  eventually  leads  to  the  formation  of  a  nonlinear 
bubble-and-spike  structure2. 

In  laser  fusion,  the  cold,  dense  pellet  is  imploded  by  the  rapid  ablation  of  hot.  low  den¬ 
sity  plasma  from  the  pellet’s  outer  surface.  This  ablative  acceleration  process  is  Rayleigh- 
Taylor  unstable.  This  desymmetrizing  instability  may  grow  at  the  ablation  layer  on  the 
outer  surface  of  the  pellet,  at  an  internal  interface  where  materials  of  different  density  abut 
or  at  the  fuel-pusher  interface  during  the  final  deceleration  phase  of  the  implosion  process: 
i.e.,  whenever  the  density  and  pressure  gradients  have  opposite  signs.  The  RT  instability 
is  a  potential  obstacle  to  laser  fusion  since  it  can  cause  shell  fracture,  fuel  mix,  or,  in  its 
mildest  form,  a  nonuniform  implosion  that  will  reduce  the  convergence  and  severly  dimin¬ 
ish  the  energy  gain  of  the  fusion  pellet.  It  is  the  ablative  RT  instability  which  we  address 
in  this  report. 

Classical,  linear  theory  predicts  that  an  initial  perturbation,  rj0,  on  the  interface  of  an 
unstably  stratified  fluid  in  planar  geometry  will  grow  as  77(f)  —  rj0eyt,  where  7  =  (Akg)1^2 
with  k  =  27r  /A.  Here  A  is  the  wavelength  of  the  perturbation,  g  is  the  acceleration.  A  = 
( pu  ~  Pi)/{Pu  +  Pi)  is  the  Atwood  number  and  pu(pi)  is  the  density  of  the  upper  (lower) 
fluid.  (A  ~  1  in  the  laser-ablative  case.)  Recent  theoretical3-10,  numerical11-17  and 
experimental18-20  investigations  of  the  RT  instability  in  an  ablative  environment  indicate 
that  the  growth  behavior  is  quite  different  from  the  predictions  of  simple  theory.  Results 
from  these  investigations  are  in  general  agreement  and  indicate  that  the  RT  growth  rate  is 
about  1/2  of  the  classically  predicted  value  for  long  wavelength  perturbations  with  1  pm 
laser  light  and  a  strong  cutoff  occurs  in  the  growth  rate  for  short  wavelength  perturbations. 

This  reduction  in  the  growth  rate  is  important  because  it  is  the  RT  instability  that 
ultimately  dictates  the  allowed  aspect  ratio  R/AR  (ratio  of  the  distance  pushed  to  the  shell 
thickness)  of  the  shell.  The  allowed  aspect  ratio  both  fixes  the  pressure  required  to  drive 
the  shell  inward  and  indirectly  determines  the  efficiency  of  the  energy  coupling  between 
the  laser  and  the  shell21.  If  the  RT  instability  grows  at  nearly  its  classically  predicted 
rate  then  thick,  low-aspect-ratio  shells  must  be  used.  Such  thick  shells  are  predominantly 
undesirable  for  two  reasons:  they  foster  the  growth  of  plasma  instabilities,  which  scale 
as  IA2  (where  I  is  the  laser  intensity  and  A;  is  the  laser  wavelength),  and  they  lead  to 
unacceptable  low  rocket  efficiencies,  in  the  range  of  2%  —  8%.  with  consequent  relatively 
low  pellet  gains  and  high  required  driver  energies  for  ignition12. 
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The  allowed  aspect  ratio  is  a  function  of  the  RT  growth  rate  and  the  wavelength  of  the 
most  dangerous  mode.  In  spherical  coordinates,  the  planar  classical  growth  rate  (kg)1/2 
goes  over  to  (liZ/R)1/2,  for  l  1.  Let  us  assume  that  the  growth  rate  for  the  laser-ablative 
case  is  e  (/iZ/R)1/2  where  e  <  1.  If  the  shell  is  imploded  halfway  inward  at  constant 
acceleration  then  the  number  of  e-foldings  of  unstable  growth22  is  -ft  =  (e7r/4)(2/)1/2. 
Classical,  incompressible  theory  suggests  that  the  most  dangerous  mode  is  A  =  A/Z,  the 
thickness  of  the  shell;  for  this  mode,  the  inner  and  outer  surfaces  are  strongly  coupled. 
However,  numerical  simulations  have  shown  that  when  compressibility  and  nonlinear  effects 
are  taken  into  account,  it  is  those  perturbation  wavelengths  approximately  3  times  the  shell 
thickness  that  are  most  damaging  to  the  implosion  process16.  For  convenience,  we  set  A 
=  2ttRq/1  as  7rAiZ,  and  the  limiting  aspect  ratio  is  given  by  iZ0/AiZ  =  (2-/t/xe)2. 

This  limit  can  be  quite  severe  since  an  initial  aspect  ratio  of  ten  or  so  can  become  an 
in-flight  aspect  ratio  (IFAR)  of  ~  150  as  a  result  of  compression.  IFAR  is  defined  as  the 
inner  radius  of  the  target  divided  by  the  distance  between  the  two  radii  where  the  density 
drops  to  1/e  of  the  peak  density.  Consider  a  pellet  polished  so  that  after  compression  the 
largest  surface  perturbation  has  an  amplitude  of  ~  50  A,  and  assume  that  the  implosion 
process  can  survive  about  seven  e -foldings  of  RT  growth.  This  gives  a  final  perturbation 
amplitude  of  ~  5.5 ^m,  which  would  be  ~  1/4  of  the  in-flight  target  thickness.  If  the 
RT  instability  grows  at  its  classically  predicted  rate  (e  =  1),  then  the  maximum  IFAR 
is  limited  to  a  value  less  than  20.  The  initial  aspect  ratio  for  this  pellet  would  be  about 
1.5,  and  the  laser  intensity  required  to  generate  sufficient  pressure  to  accelerate  this  shell 
to  implosion  velocities  would  be  well  above  the  threshold  for  plasma  instabilities.  These 
thick  low-aspect-ratio  shells  would  also  severely  limit  the  efficiency  and  gain  of  the  fusion 
pellet21. 

As  noted  above,  there  is  general  agreement  that  the  RT  growth  rate  is  about  1/2  of 
the  classically  predicted  value  (c  =  1/2)  with  1/im  laser  light.  This  would  permit  an  IFAR 
~  77,  which  is  marginal  at  best.  Recent  results9,1 1  have  indicated  that  the  RT  growth 
rate  is  reduced  to  about  30%  of  the  classical  value  with  1/4 /im  laser  light.  With  e  =  0.30. 
the  maximum  IFAR  «  220.  These  results  imply  that  targets  with  an  IFAR  ~  200  can 
be  stably  imploded.  Both  the  hydrodynamic  efficiency  and  maximum  laser  intensity  are 
strong  functions  of  the  IFAR22.  With  an  IFAR  of  200.  efficiencies  of  greater  than  15% 
are  attainable  with  laser  intensities  of  a  few  times  1014W/cm2.  This  further  implies  pellet 
gains  greater  than  200,  which  suggests  that  high-aspect-ratio  shells  are  a  viable  design  for 
direct-drive  laser  fusion. 

It  is  the  strong  dependence  of  the  RT  growth  rate  on  the  laser  wavelength  which 
we  address  in  this  paper.  Since  the  above  analysis  assumes  perturbation  wavelengths 
small  compared  to  the  target  radius  (large  /  numbers),  we  investigate  the  evolution  of 
the  RT  instability  in  planar,  laser  ablatively  accelerated  targets.  We  present  results  from 
our  FAST2D  Laser-Matter-Interaction  Code  for  a  series  of  simulations  which  address  the 
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issues  of  target  thickness,  perturbation  wavelength,  laser  intensity  and  laser  wavelength. 
The  results  are  compared  to  our  theoretical  model  where  the  inhibited  growth  is  based 
on  the  ablative  convection  of  vorticity  away  from  the  unstable  ablation  layer.  In  the  next 
section  we  describe  the  numerical  model  we  use  to  perform  the  simulations  and  discuss  the 
initial  conditions. 

II.  MODEL 

The  evolution  of  the  Rayleigh-Taylor  instability  in  laser  ablatively  accelerated  tar¬ 
gets  is  modeled  with  our  FAST2D  Laser  Matter  Interaction  Code.  This  is  a  fully  two- 
dimensional,  Cartesian,  hydrodynamics  code  with  a  sliding  Eulerian  grid  with  variable 
grid  spacing.  The  grid  spacing  is  typically  0.25  nm  (  0.125  pim  for  the  high  laser  intensity 
cases  )  for  15  -  20  zones  on  either  side  of  the  ablation  layer  and  increases  uniformly  to  2  - 
10  pm  for  most  of  the  rest  of  the  grid.  The  grids  at  the  outer  boundaries  of  the  compu¬ 
tational  mesh  are  stretched  to  200  -  300  /jm.  The  grid  slides  with  the  average  velocity  of 
the  accelerated  foil  so  that  the  ablation  front  remains  approximately  centered  in  the  finely 
zoned  region  thus  reducing  convective  errors. 

FAST2D  is  fully  Eulerian  which  is  necessitated  by  the  requirement  to  accurately  model 
rotational  flows.  It  should  be  pointed  out  that  the  flows  stemming  from  the  RT  instability 
are  inherently  rotational.  Vorticity  is  generated  at  the  unstable  ablation  front  by  the 
baroclinic  (noncollinear  density  and  pressure  gradients)  nature  of  the  profiles.  The  vorticity 
that  arises  from  the  differential  twisting  of  the  density  and  pressure  surfaces  feeds  the 
“heavy”  material  into  the  spike  and  the  “light”  material  into  the  bubble.  Rectilinear 
Lagrangian  algorithms  cannot  model  rotational  flows  accurately  without  restructuring  the 
grid  which  typically  introduces  an  unacceptable  amount  of  numerical  diffusion. 

The  time-dependent  equations  for  the  continuity  of  mass,  momentum  and  total  energy 
sure  integrated  numerically  with  sliding-zone  flux- corrected-  transport  (FCT)23.  FCT  solves 
the  continuity  equations  for  the  fluxes  relative  to  moving  zone  boundaries.  The  overall 
scheme  is  second-order  accurate  and  no  artificial  viscosity  is  explicitly  required  to  maintain 
numerical  stability.  Mass,  momentum  and  energy  are  both  locally  and  globally  conserved. 

Thermal  conduction  is  incorporated  in  a  single- temperature  model  using  a  semi- 
implicit  solution  of  the  two-dimensional,  classical,  non-linear  thermal  conduction  equation 
based  on  the  model  of  Spitzer  and  Harm24.  The  laser  energy  is  absorbed  by  classical  in¬ 
verse  bremsstrahlung  absorption  with  a  resonant  dump  at  the  critical  density.  The  laser 
light  is  at  normal  incidence.  The  fluid  equations  are  closed  with  an  ideal  equation  of  state. 

FAST2D  and  its  counterpart  in  cylindrical  geometry,  FASTRZ.  have  been  extensively 
compared  with  experimental  results  and  have  been  found  to  give  accurate  results  over  a 
wide  range  of  parameter  space.  These  comparisons  include  scaling  laws  (mass  ablation 
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rates,  target  velocities,  velocities  of  ablated  material  and  ablation  pressures)25,26,  velocity 
profiles  of  the  ablated  material27,  target  velocity  asymmetries  driven  by  imposed  beam 
nonuniformities28,  rear  surface  temperatures  due  to  shock  heating  and  thermal  transport29, 
laser  absorption  rates  and  finite  spotsize  effects30,  temperature  and  density  profiles  in 
the  underdense  ablated  plasma31  and  the  linear20  and  nonlinear19  evolution  of  the  RT 
instability  in  structured  targets. 

In  order  to  obtain  well  defined  RT  growth  rates  we  utilize  a  new  “quiet  start”  method. 
FAST2D  is  run  in  a  one-dimensional  mode  until  a  nearly  steady-state  distribution  is  at¬ 
tained.  We  start  with  a  laser  pulse  which  has  a  2  ns  HWHM  (half-width-half-maximum) 
Gaussian  rise  after  which  the  laser  intensity  is  held  constant  at  the  peak  value.  A  steady- 
state  is  typically  reached  ~  2  ns  after  the  peak  of  the  laser  pulse.  These  steady-state 
profiles  then  serve  as  initial  conditions  in  the  two-dimensional  model.  At  the  start  of  each 
RT  simulation,  these  steady-state  density  profiles  axe  perturbed  at  the  ablation  front  with 
a  single  wavelength  perturbation.  Small  sinusoidal  variations  in  the  density  are  impressed 
in  10  zones  on  the  high  density  side  of  the  ablation  layer.  Initial  perturbation  strengths 
vary  slightly  and  are  bounded  from  above  by  a  total  mass  perturbation  of  less  than  0.5%. 
The  growth  rate  of  the  instability  is  then  measured  during  the  constant  intensity  phase  of 
the  laser  pulse  during  which  time  the  target  acceleration  is  approximately  constant.  This 
technique  has  several  advantages:  it  avoids  the  mixing  of  RT  phenomena  and  simple  sound 
waves  bouncing  off  the  density  perturbations;  we  are  assured  of  having  a  self-consistent 
steady  state  profile  which  is  difficult  to  obtain  analytically  with  inverse  bremsstrahlung 
absorption;  and,  we  save  computer  time  by  using  the  same  initial  conditions  for  each 
perturbation  wavelength,  for  fixed  laser  wavelength  and  intensity. 

The  computational  growth  rates  are  obtained  by  Fourier  transformation  of  the 
summed  mass  of  the  foil.  The  mass  of  the  foil  {pR)  is  integrated  from  the  rear  of  the 
foil  (the  side  away  from  the  laser)  to  the  ablation  front  for  each  transverse  coordinate.  A 
plot  of  the  time  evolution  of  the  Fourier  transform  coefficients  provides  a  very  sensitive 
measure  of  the  increase  in  areal  mass  nonuniformities  due  to  the  RT  instability.  Typically 
5-7  e-foldings  of  exponential  growth  are  obtained  for  the  results  presented  here. 

The  foils  used  in  this  study  are  plastic  (CH).  The  initial  foil  thicknesses  are  55  pm 
and  150  pm  with  absorbed  laser  intensities  of  7  x  1013  W/cm2  and  3  x  1014  W/cm’ 
respectively.  The  evolution  of  the  RT  instability  on  both  targets  is  investigated  with  laser 
wavelengths  of  1,  1/2  and  1/4  pm.  Perturbation  wavelengths  range  from  5  pm  to  150  pm. 
The  numerical  results  are  presented  in  the  next  section. 

III.  NUMERICAL  RESULTS 

We  first  compare  the  evolution  of  the  RT  instability  with  and  without  the  “quiet  start" 
technique.  The  results  are  shown  in  Fig.  1.  Here  is  plotted  the  loge  of  the  amplitude  of 


the  Fourier  transform  coefficient  for  the  perturbed  mode  (A  =  2 O^im)  as  a  function  of 
time  for  both  cases.  The  initial  target  thickness  is  55  nm  with  1  /um  laser  light  at  7  x 
1013  W/cm2.  The  square  symbols  are  for  the  case  where,  starting  from  t  =  0,  the  laser 
has  a  2  ns  Gaussian  rise  after  which  the  intensity  is  held  constant.  There  is  an  initial 
transition  phase  of  approximately  0.75  ns  duration  while  the  flow  patterns  evolve  into  an 
eigenmode  followed  by  rapid  growth.  After  several  oscillations,  the  growth  finally  evolves 
into  a  quasi-uniform  exponential  growth  phase.  The  initial  rapid  oscillatory  growth  is  most 
likely  due  to  a  combination  of  the  strong  shocks  generated  by  the  rising  laser  pulse  and  a 
non-self-consistent  initial  steady-state.  The  triangles  axe  for  the  “quiet  start”  case  and  the 
t=0  point  is  actually  4  ns  after  the  laser  is  "turned  on”.  The  laser  intensity  is  constant 
throughout  the  time  plotted;  and,  after  a  slightly  shorter  transition  period  (0.50  ns)  the 
growth  is  quite  uniform.  The  growth  rate  is  1.49  x  109  ns-1.  The  “quiet  start”  method 
is  used  for  all  the  results  presented  here. 

A  comparison  of  the  numerical  growth  rates  (square  symbols)  to  the  classical  growth 
rates  (triangle  symbols)  as  a  function  of  perturbation  wavelength  for  1  /j.m  laser  light  is 
shown  in  Fig.  2.  Fig.  2a  shows  the  results  for  55  fim  thick  CH  at  7  x  1013  W/cm2  and 
Fig.  2b  is  for  150  fxm  thick  CH  at  3  x  1014  W/cm2.  In  both  cases  the  computational 
growth  rates  follow  an  approximately  k1/2  dependence  with  a  moderately  strong  cutoff  for 
perturbation  wavelengths  less  than  10  gm.  For  perturbation  wavelengths  in  the  50  fim 
range,  the  computational  growth  rate  is  reduced  to  ~  75%  of  the  classical  growth  rate 
for  the  moderate  intensity  case  (7  x  1013  W/cm2)  and  to  ~  65%  of  classical  for  the  high 
intensity  case  (3  x  1014  W/cm2).  The  classical  growth  rate  is  given  by  (kg)1/2  where  g  = 
P„/M.  Pa  is  the  average  ablation  pressure  across  the  foil  surface  and  M  is  the  mass  density 
of  the  foil.  For  the  moderate  intensity  case  Pa  =  8.9  Mb  and  M  =  4.94  x  10~3  gm/cm2 
and  thus  g  =  1.8  x  1015  cm/s2.  For  the  high  intensity  case  Pa  =  24.5  Mb  with  M  =  1.23 
x  IQ-2  gm/cm2  and  g  =  2.0  x  1015  cm/s2.  In  each  case  there  were  only  slight  variations 
in  g,  0(±1%),  as  the  perturbation  wavelength  was  varied. 

The  results  with  1/2  fim  laser  light  are  shown  in  Fig.  3.  Again,  the  moderate  intensity 
comparison  is  in  Fig.  3a  and  the  high  intensity  comparison  is  in  Fig.  3b.  In  both  cases 
the  computational  growth  rates  are  more  strongly  inhibited  than  in  the  i  fim  case  and  the 
dependence  on  the  perturbation  wavelength  is  reduced.  The  short  perturbation  wavelength 
cutoff  is  not  as  dominant  as  in  the  1  pm  case  or  it  occurs  at  wavelengths  less  than  5  /jm. 
The  growth  rate  is  reduced  to  50  %  (45%)  of  classical  for  the  moderate  (high)  intensity 
case  for  wavelengths  0(50  f. im ).  The  parameters  are  Pa  =  10  Mb.  M  =  5.0  x  10-3  gm/cm2 
and  g  =  2.0  x  1015  cm/s2  for  the  moderate  intensity  case:  and.  Pa  =  30.5  Mb.  M  =  1.27 
x  10~2  gm/cm2  and  g  =  2.4  x  1015  cm/s2  for  the  high  intensity  case. 

The  RT  growth  rates  are  even  more  strongly  inhibited  with  1/4  fim  laser  light  as 
shown  in  Fig.  4.  The  RT  growth  rate  is  strongly  reduced  for  both  the  moderate  ( Fig.  4a) 
and  high  (Fig.  4b)  intensity  cases  and  the  dependence  on  the  perturbation  wavelength  is 


strongly  suppressed.  For  A  ~  50  pm,  the  growth  rate  is  only  30%  of  the  classical  value. 
For  the  moderate  intensity  case,  P„  =  15.5  Mb,  M  =  4.S  x  10-3  gm/cnr  and  g  =  3  2  x 
1015  cm/s2;  and,  for  the  high  intensity  case,  P„  =  40.6  Mb,  M  =  1.25  x  10~J  gm/cnr  and 
g  =  3.25  x  1015  cm/s2. 

These  results  are  summarized  are  Fig.  5  where  the  ratios  of  the  numerical  growth 
rates  to  the  classical  growth  rates  axe  plotted  as  a  function  of  the  perturbation  wavelength 
and  the  laser  wavelength  for  both  the  moderate  (Fig.  5a)  and  high  (Fig.  5b)  intensity 
laser  cases.  It  is  quite  clear  that  the  RT  growth  rate  is  strongly  dependent  upon  the  laser 
wavelength.  With  the  exception  of  the  1  pm  results,  there  appears  to  be  little  dependence 
upon  the  laser  intensity.  Resonant  absorption  is  more  dominant  with  1  pm  laser  light  and 
this  may  be  influencing  the  results.  This  is  addressed  in  more  detail  in  the  next  section. 
The  results  indicate  that  the  growth  rate  reduction  may  be  asymptoting  to  a  constant 
reduction  factor  for  the  very  long  wavelength  modes;  see  Fig.  5b.  The  reduction  factors 
are  1/1.54,  1/1.82  and  1/3  for  the  very  long  wavelength  perturbations  for  1,  1/2  and  1/4 
pm  laser  light  respectively.  In  the  next  section  the  theoretical  model  is  summarized  and 
compared  to  the  numerical  results. 

IV.  THEORETICAL  RESULTS 

There  have  been  many  theoretical  models  developed  over  the  past  several  years  in  an 
attempt  to  explain  the  inhibited  RT  growth  in  ablatively  accelerated  targets.  These  models 
have  proposed  convection3,  mass  ablation5,6,9,  thermal  conduction3,5,  compressibility4  •t>, 
combinations  of  the  above6  and  density  gradient  stabilization7  as  the  cause  of  the  inhibited 
growth.  Recently  it  has  been  shown  that  compressibility  plays  a  minor  role4,8  and  can. 
in  fact,  be  destabilizing8.  In  direct-drive  laser  fusion  the  density-gradient  scalelengths  are 
too  small  to  significantly  reduce  the  RT  growth  rate  by  density-gradient  stabilization1 ' . 
Fig.  6  compares  the  steady-state  density  and  pressure  profiles  for  both  the  1  pm.  at  10 13 
W/cm2,  and  1/4  pm,  at  3  x  1014  W/cm2,  cases.  The  density-gradient  scalelengths  are 
0.49  and  0.40  pm  respectively  and  in  neither  case  would  the  density  gradients  provide  a 
significant  reduction  in  the  growth  rate. 

The  vorticity  generated  by  the  noncollinear  nature  of  the  density  and  pressure  profiles 
is  the  controlling  element  in  the  evolution  of  the  RT  instability1, .  In  particular  it  is  the 
ablative  convection  of  vorticity  away  from  the  unstable  ablation  layer  that  reduces  the  rate 
at  which  the  fluids  are  interchanged.  The  details  of  the  analytic  scheme  will  be  published 
elsewhere32  and  we  will  only  summarize  the  results  here. 

We  consider  the  oscillations  of  an  inviscid.  thin  layer  of  fluid  (of  thickness,  h)  with  an 
exponential  density  gradient  bounded  on  both  sides  by  inviscid  fluids  of  constant  density 
of  semi-infinite  extent.  In  this  two-dimensional  geometry,  the  equation  describing  the  time 
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evolution  of  the  vorticity  of  the  thin  layer  is  obtained  by  taking  the  curl  of  the  equation 
of  motion  and  is  given  by 


Q  =  Vp  x  Vp/p2  —  V  ■  VuJ  —  lJ( V  ■  V), 


av  - 1) 


where  p  is  the  density  and  p  is  the  pressure.  The  first  term  on  the  RHS  is  the  baroclinic 
source  term  which  will  generate  vorticity  whenever  the  density  and  pressure  gradients  are 
noncollinear.  See  Fig.  6.  The  region  between  the  density  peak  and  the  pressure  peak  is 
the  RT  unstable  region33.  The  second  and  third  terms  axe  the  convective  and  compressible 
terms  respectively. 

We  first  consider  the  incompressible,  non-ablative  case  and  perform  a  first  order  per¬ 
turbation  expansion  on  the  above  equation.  Expressing  the  perturbed  velocities  in  terms 
of  the  stream  function  and  evaluating  the  boundary  conditions,  the  growth  rate  for  the 
interfacial  mode  takes  the  form 


2„  [AL  +  kh{AL  +  2)]„  ,  r,  ikh(L  +  Akh)(A  +  L),'/2 
7A9=  +  + 


(IV  —  2) 


where  A  =  (pu  -  pf)/(p„  +  pi)  is  the  Atwood  number  with  L  =  In (pu/pi)  and  pu(pi)  is  the 
density  of  the  upper  (lower)  fluid,  ptt  =  pi  eK!H .  Eq.  (IV-2)  can  be  simplified  by  setting 
L  =  h/H,  and  for  kh  small 


7 2 /kg  .4(1  +  kh)(l  —  kh/A2),or 


I2 /kg  as  1/(1  +  k2h2),  (IV-  3) 

for  A  — *  1.  The  mode  is  unstable  if  p„  >  p/  and  the  intermediate  layer  reduces  the  growth 
rate  from  the  classical  value.  The  classical  value  is  recovered  in  the  limit  kh  — *  0.  In 
the  typical  laser  ablativelv  accelerated  target,  the  thickness  of  the  ablation  region  is  quite 
small,  h  ~  0(1-2  pm)  (see  Fig.  6i,  and  there  will  be  little  reduction  in  the  growth  rate 
except  for  very  short  wavelength  perturbations. 

To  account  for  the  ablative  convection  of  vorticity  away  from  the  unstable  layer,  we 
incorporate  an  exponentially  varying,  zeroth-order  ablation  velocity  (  V’a  )  which  satisfies  the 
steady-state  continuity  equation.  V  ■  i  pi’,,  i  =  0.  Assuming  incompressible  perturbations. 


7 


we  obtain  a  fifth  order  dispersion  relation.  In  the  limit  kh  — >  0.  this  reduces  to  a  quadratic 
and  the  growth  rate  for  the  infinitesimally  thin  layer,  ablative  RT  case  is 


IV  - 


where  70  =  [Akg)1^2  is  the  classical  result.  This  expression  is  similar  those  in  Ref.  3  and  C 
except  that  Ref.  3  uses  the  target  velocity  (and  thus  little  stabilization)  and  Ref.  G  invokes 

supersonic  flow. 

If  the  ablation  velocity  is  large.  Eq.  IV-4  predicts  strongly  inhibited  growth  of  the 
instability;  the  difficulty  is  in  what  to  choose  for  the  ablation  velocity.  Figure  7  shows 
the  isovorticity  contours  on  the  target  surface  early  in  the  evolution  of  the  20  pm  RT 
instability  for  1/4 pm,  7  x  1013  W/cm2  laser  light  on  a  55  pm  thick  plastic  target.  The 
laser  light  is  impinging  on  the  target  from  the  right.  The  ablation  velocity  (V0.  cm/si, 
which  is  convecting  vorticity  away  from  the  target  surface,  is  plotted  on  top  of  the  vorticitv 
contours.  In  order  to  compare  the  above  theoretical  expression  with  the  numerical  growth 
rates,  the  ablation  velocities  are  taken  from  the  numerical  simulations  and  are  measured 
at  the  center  of  the  unstable  shear  layer  on  the  ablation  front. 

The  theoretical  growth  rates  (solid  circles)  are  compared  to  the  numerical  growth  rates 
in  Figures  2-4.  The  agreement  is  quite  good  for  the  long  wavelength  modes.  The  lack  of 
agreement  at  the  short  wavelength  perturbations,  especially  for  1  /im  laser  light,  is  due  to 
the  assumption  of  an  infinitely  thin  vortex  sheet  at  the  ablation  front.  Eq.  IV-4  cannot 
reproduce  the  short  wavelength  cutoff.  In  order  to  account  for  the  finite  size  nature  of  the 
unstable  layer,  we  replace  the  classical  expression  in  Eq.  IV-4  (70)  by  the  approximate 
expression  for  the  growth  rate  for  an  exponential  density  layer  of  finite  thickness  (Eq.  IV- 
3).  These  results  are  shown  as  the  open  circular  symbols  in  Figs.  2-4.  This  modification 
to  Eq.  IV-4  has  little  effect  on  the  long  wavelength  modes  and  does  a  reasonable  job  of 
reproducing  the  short  wavelength  cutoff,  h  is  taken  to  be  the  distance  between  the  density 
peak  and  the  pressure  peak  as  determined  from  the  steady-state  profiles.  The  size  of  h  is 
laser  intensity  dependent;  h  =  1.0  pm  for  the  moderate  intensity  case  and  0.75  pm  for  the 
high  intensity  case. 

The  strong  reduction  in  RT  growth  with  short  wavelength  lasers  is  due  to  the  increase 
in  ablation  velocity  as  the  laser  wavelength  is  reduced.  The  higher  mass  ablation  rate  for 
short  wavelength  lasers  produces  a  higher  acceleration  of  plasma  away  from  the  ablation 
surface  and  therefore  a  larger  value  of  Va.  Figure  8  is  a  plot  of  the  ablation  velocity  i  V,t  1 
obtained  from  the  FAST2D  Model  as  a  function  of  distance  from  the  target  surface  for  1. 
1/2  and  1/4  pm  laser  light.  All  the  calculations  were  with  a  55  pm  thick  CH  target  with 
an  absorbed  intensity  of  7  x  1013  W/cm2  under  nearly  steady-state  conditions.  Note  that 
at  the  center  of  the  unstable  shear  layer,  the  ablation  velocity  for  1/4  pm  light  is  over 


three  times  larger  than  for  1  pm  light  These  re-::l*s  are  -.ub>*ant:ared  by  simple  scaling 
arguments* '  which  show  that  Y„  >  \7l  The  growth  rate  for  the  RT  instability  will  be 
more  strongly  inhibited  with  short  wavelength  lasers  than  w>th  1  pm  laser  light,  since  the 
vorticity  will  be  converted  away  at  a  faster  rate  Uf  course,  at  large  distances  from  the 
ablation  layer  the  blowoff  velocity  is  less  for  shorter  wavelength  lasers 

It  is  somewhat  difficult  to  obtain  a  definitive  scaling  relationship  between  the  ablation 
velocity  and  the  laser  wavelength  and  intensity  This  is  because  of  the  dependence  of  the 
thickness  of  the  unstable  layer  on  the  laser  intensity  and  the  strongly  exponential  nature 
of  tiie  ablation  velocity.  For  low  intensity  lasers.  .  10n  \Y  cm*  i.  the  distance  between  the 
density  peak  and  pressure  peak  is  --  1.5  pm  and  the  center  of  the  vortex  layer  is  nearly 

2  pm  from  the  density  peak  •  see  Ref.  16  As  a  result,  at  this  low  intensity,  the  ablation 
velocity  through  the  vortex  sheet  is  quite  large.  >  LO*  cm/s.  For  the  moderate  and  high 
intensity  cases  simulated  here,  the  distance  between  the  density  and  pressure  peaks  is  1 
pm  or  less,  and  the  center  of  the  unstable  layer  is  about  a  micron  from  the  density  peas. 
This  results  in  an  ablation  velocity  through  the  vortex  layer  of  less  than  10*  cm/s  for  1  pm 
light  at  7  x  10i3  W  cm*  Note  however,  at  comparable  distances  1  ~~  2  pm  i,  the  ablation 
velocity  is  greater  than  that  for  1013  W  cm*  Further  work  is  in  progress  to  investigate 
this  ablation  layer  scaling. 

Recent  work  by  Takabe.  et  al.9  also  showed  an  enhanced  reduction  in  the  RT  growth 
rate  with  1/4  pm  laser  light.  There  they  applied  am  eigenvalue  analysis  to  the  linearized 
equations  for  stationary  ablating  plasmas.  Their  results  agreed  well  with  ours13  and  those 
of  McCrorv  et  al.12  for  1  pm  laser  light.  The  growth  rates  they  obtained  with  1/4  pm 
light  were  factors  of  2-3  smaller  than  the  rates  obtained  with  the  FAST2D  Model.  The 
most  likely  source  of  this  discrepancy  is  due  to  their  assumption  that  all  of  the  laser  light 
is  absorbed  at  the  critical  density.  We  find  in  the  numerical  simulation  that  most  of  the 
laser  energy  is  absorbed  through  inverse  bremsstrahlung  with  1/4  pm  light  and  thus  the 
temperature  profile  is  not  nearly  so  steep  as  that  resulting  from  resonant  absorption.  As 
a  result,  the  mass  ablation  rate  with  inverse  bremsstrahlung  absorption  is  not  as  large  as 
one  would  expect  from  scaling  laws  based  or.  energy  deposition  at  critical  density  and  the 
growth  rate  is  not  reduced  by  as  large  a  factor. 

To  test  this  hypothesis,  we  ran  FAST2D  with  100VT  deposition  at  critical  for  1/4  pm. 

3  x  1014  W /cm1  laser  light  on  a  150  pm  thich  CH  target.  Under  these  conditions,  the 
ratio  of  numerical  to  classical  growth  rate  is  0.2  which  is  2/3  of  the  growth  rate  obtained 
with  inverse  bremsstrahlung  absorption.  The  analytical  expression  obtained  by  Takabe  et 
al.  is  based  on  the  initial  target  mass.  If  this  is  modified  fo  accommodate  the  accelerated 
target  mass,  which  is  ~  50%  of  the  initial  target  mass,  their  growth  rates  would  increase 
bv  10  -15  %  and  would  be  in  good  agreement  with  our  resonant  absorption  results.  We 
believe  that  these  two  effects  explain  the  very  strong  reduction  in  growth  rates  with  1,4 
pm  laser  light  obtained  by  Takabe.  et  al.9 
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V.  NONLINEAR  RESULTS 


As  corrugations  on  the  interface  of  an  unstably  stratified  fluid  grow  in  amplitude,  they 
eventually  form  the  nonlinear  bubble- and-spike  structure2.  The  spikes  are  characterized  by 
thin  sheets  (in  2  dimensions)  of  dense  fluid  which  “fall”  and  the  bubbles  are  rounded  regions 
of  light  fluid  which  “rise”.  When  the  amplitude  of  the  spike  becomes  large  ,  0(A/2),  the 
spike  partially  shields  the  bubble  from  the  ablation  process  and  the  vorticity  is  advected 
down  the  sides  of  the  spike  causing  a  Kelvin- Helmholtz-like  (KH)  rollup  of  the  spike  tips. 
We  have  previously  reported  on  this  phenomena  with  low  intensity  laser  light12’15.  We 
find  a  somewhat  similar  process  occuring  with  high  intensity  lasers. 

Figure  9a  shows  the  isodensity  contours  after  ~  4  e -foldings  of  growth  for  1/4  pm 
light  at  7  x  1013  W/cm2.  The  wavelength  of  the  perturbation  is  20  pm.  At  this  time  the 
spike  region  is  still  quite  broad  and  the  bubble  region  is  relatively  flat.  The  corresponding 
streamlines  are  plotted  in  Figure  9b.  Note  the  strong  vortical  flow  at  the  base  of  the  spike. 
The  flow  pattern  (in  a  case  for  which  all  the  fluid  flow  is  from  left  to  right)  is  quite  similar 
to  the  vortex  ring  that  is  formed  when  fluid  eminates  from  a  nozzle  .  Figure  10  shows  the 
isodensity  contours  at  a  slightly  later  time.  The  bubble  is  now  rounded  and  the  spike  is 
much  narrower  with  significant  broadening  of  the  spike  tip. 

Figures  11  and  12  illustrate  the  nonlinear  bubble-and-spike  structure  for  1/2  pm,  7 
x  1013  and  1  pm,  3  x  1014  laser  light  respectively.  Although  all  the  above  cases  indicate 
some  degree  of  KH-like  rollup  of  the  spike  tips,  the  degree  of  rollup  is  not  as  large  as  that 
obtained  for  low  intensity  laser  light  and  the  bubble  region  remains  quite  flat.  This  is 
due  to  the  higher  mass  ablation  rate  with  the  higher  intensity  and/or  shorter  wavelength 
lasers.  This  is  more  in  agreement  with  the  work  of  McCrory  et  al.12  find  Verdon  et  al.14 
than  was  our  previous  lower  intensity  case13.  Note  that  even  with  the  relatively  large  spike 
amplitude,  the  bulk  of  the  target  is  still  flat  and  is  uniformly  accelerated. 

VI.  CONCLUSIONS 

We  have  shown,  through  a  series  of  detailed  numerical  simulations,  that  the  Ravleigh- 
Tavlor  growth  rate  for  laser  ablatively  accelerated  targets  is  a  strong  function  of  laser 
wavelength.  The  numerical  results  agree  well  with  the  theoretical  model  based  on  the 
convection  of  vorticity  away  from  the  unstable  ablation  front.  As  the  laser  wavelength  is 
reduced,  the  ablation  velocity  increases  thereby  convecting  vorticity  away  at  a  faster  rate 
and  reducing  the  rate  at  which  the  fluids  are  interchanged. 

If  we  consider  perturbation  wavelengths  of  ~  60  pm.  the  RT  growth  rate  is  reduced 
by  factors  of  1/1.54,  1/2.22  rind  1/3.33  for  1,  1/2  and  1/4  pm  laser  light  respectively. 
Assuming  the  shell  can  survive  up  to  7  e-foldings  of  exponential  growth  at  this  perturbation 
wavelength,  the  maximum  in-flight  aspect  ratio  for  a  reactor  size  pellet  is  R/AR  47. 


98  and  220  for  1,  1/2  and  1/4  fim  laser  light.  It  has  been  shown22  that  with  an  IFAR  of 
200,  efficiencies  of  greater  than  15%  are  attainable  with  1/4  ytm  light  and  intensities  of 
a  few  times  1014  W/cm2.  The  results  presented  here  imply  that  targets  with  an  initial 
aspect  ratio  of  ~  10  and  an  IFAR  ~  200  can  be  stably  imploded.  This  suggests  that 
moderately  high- aspect-ratio  shells  are  a  viable  design  for  direct-drive  laser  fusion  with 
short  wavelength  (1/4  ^m)  laser  light. 

It  is  clear  that  ablation  layer  phenomena  dominates  the  development  and  evolution 
of  the  RT  instability  in  laser  ablatively  accelerated  targets.  The  convection  velocity  at 
the  ablation  front  controls  the  degree  of  inhibition  of  the  RT  growth.  The  thickness  of 
the  ablation  front  strongly  influences  the  cutoff  at  short  perturbation  wavelengths.  More 
work  is  needed  in  this  area,  to  develop  scaling  relationships  for  ablation  layer  parameters 
as  a  function  of  laser  wavelength  and  intensity  and  target  composition,  all  under  inverse 
bremsstrahlung  absorption  conditions.  We  also  need  additional  experimental  data  on  the 
RT  instability  to  test  these  predictions. 
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Fig.l.  Plot  of  the  log*  of  the  Fourier  transform  coefficient  of  the  Rayleigh- Taylor  unstable 
mode  as  a  function  of  time  for  both  the  nonsteady-state  (squares,  temporally  rising  laser 
pulse)  and  steady-state  (triangles,  “quiet  start”,  constant  laser  pulse)  initial  conditions. 
In  both  cases  the  initial  target  thickness  was  55  pm  with  1  pm  laser  light  with  a  peak 
intensity  of  7  x  1013  W/cm2  and  a  perturbation  wavelength  of  20  pm. 
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Fig. 5.  The  ratios  of  the  numerical  growth  rates  to  the  classical  growth  rates  as  a  function 
of  perturbation  wavelength  and  laser  wavelength.  The  triangle,  circle  and  square  symbols 
are  for  1,  1/2  and  1/4  fim  laser  light  respectively.  Fig.  5a  is  for  an  intensity  of  7  x  1013 
W/cm2  and  Fig  5b  is  for  an  intensity  of  3  x  1014  W/cm2. 


Fig. 7.  Isovorticity  contours,  in  12%  increments  of  the  maximum  on  the  target  surface  early 
in  the  evolution  of  the  20  pm  RT  instability  with  1/4  pm,  7  x  1013  W/cm:  laser  light.  The 
laser  is  impinging  the  target  from  the  right.  The  spatial  dimensions  axe  in  microns.  The 
dashed  (solid)  contours  indicate  clockwise  (counterclockwise)  flow.  The  ablation  velocity 
(V0,  cm/s)  is  plotted  on  top  of  the  vorticity  contours.  The  solid  vertical  line  illustrates 
the  point  at  which  the  ablation  velocity  is  measured. 
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